Helical modes in carbon nanotubes generated by strong electric fields 
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Helical modes, conducting opposite spins in opposite directions, are shown to exist in metallic 
armchair nanotubes in an all-electric setup. This is a consequence of the interplay between spin- 
orbit interaction and strong electric fields. The helical regime can also be obtained in chiral metallic 
nanotubes by applying an additional magnetic field. In particular, it is possible to obtain helical 
modes at one of the two Dirac points only, while the other one remains gapped. Starting from a 
tight-binding model we derive the effective low-energy Hamiltonian and the resulting spectrum. 
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Carbon based solid state physics has attracted much 
attention over the past decades. One of the best studied 
structures in this field is the carbon nanotube (CNT), a 
hexagonal lattice of carbon atoms rolled up to a cylinder 
[l|. The experimental techniques for creating, isolating, 
and analyzing CNTs have by now remarkably matured, 
such that characteristics that have previously been ob- 
scured by disorder can now be experimentally resolved 
[2-5]. An example is the spin-orbit interaction (SOI), 
which is generally small in CNTs yet can affect 

electron spin decoherence in CNT quantum dots |10L [ll| , 
or allow spin control 



12 



H3| and spin filtering pj. A 
complete understanding of the SOI in CNTs becomes 
therefore desirable. 

In this Letter, we investigate the effect of SOI in com- 
bination with a strong electric field in single-wall CNTs 
within an effective low-energy theory. In particular, we 
identify experimentally accessible parameter regimes in 
which SOI and electric fields create helical modes with- 
out the need for magnetic fields. This must be contrasted 
with the helical modes in one-dimensional metals with 
Rashba SOI, which can be created only with an addi- 
tional magnetic field that opens a gap at the crossing 
point of the two Rashba-shifted parabolas 14, lB| . He- 
lical modes, conduction channels transporting opposite 
spins in opposite directions, naturally lead to spin filter- 
ing, but they have also potential application as Cooper 
pair splitters [16| and, if in proximity with a supercon- 
ductor, lead to Majorana bound states at the edges of 
the conductor [ItJ • Helical modes have also attracted 
much attention recently in the context of topological in- 
sulators [18[. Such physics may be achieved in CNTs in 
an all-electric setup. 

Perfect helical modes appear in armchair CNTs, while 
in metallic chiral CNTs the spins of the left and right 
moving modes are not precisely opposite. In the latter, 
however, perfect helicity can be restored in one Dirac 
point by an additional magnetic field, whereas the other 
Dirac point becomes insulating at these energies. This 
corresponds to the effective suppression of one valley for 
the low-energy physics. 

The model The effective theory is based on a compre- 




FIG. 1. (Color online) Cross section of a CNT in a uniform 
electric field E. The orientation of the orbitals p r ,pt as well 
as the local coordinate system r, t depends on the azimuthal 
angle ip. The s orbital is indicated by the dashed circles. The 
electric field E is oriented along the x-direction of the global 
coordinate system. The z-direction is along the nanotube. 

hensive model which incorporates the curvature effects 
for nearest-neighbor hopping and orthogonal orbitals [l[ . 
Charge effects in CNTs due to electric fields have been 
considered before [l9l-[2lj. Here, we also include spin ef- 
fects induced by external uniform electric fields (see Fig. 
[I]). For this we start from a tight-binding description of 
the honeycomb lattice on a cylinder surface where we in- 
clude all orbitals of the second shell and the hybridization 
of the 7r and the a bands. The screening of the electric 
field by electron-electron interactions is treated on the 
mean-field level. The corresponding Hamiltonian is 



H = H hs + Hso + H 



(i) 



H 



(2) 



(1) 



The band structure Hamiltonian H\y S includes the hop- 
ping of electrons between orbitals of neighboring car- 
bon atoms and accounts for the orbital energies = 

^tf c ifjL\ c 3t*'* + £sc\ sX c is \. Here c ifI \ are the electron op- 
erators, i and j are nearest neighbor sites on the hon- 
eycomb lattice, A = ±1 is the spin in z-direction, and 
\i runs over the second shell orbitals with \i = s the s 
orbital and \i = p r -,Pt-,Pz the p orbitals pointing in ra- 
dial, tangential and ^-direction (see Fig. [I]). The tt band 
is formed by the p r orbitals, while the a band is formed 
by PtiPz-, s - Summation over repeated indices is assumed. 
The hopping amplitude between (j, //) and (i, \i) is a 
linear combination of the four fundamental hopping am- 
plitudes V ss , V sp: Vp p , Vp p [1] with coefficients depending 
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on the relative orientation of the orbitals ji and \j! [22j. 
The energy difference between s and p orbitals is e s . 
The atomic SOI is modeled by the on-site Hamilto- 

where now //, rj = 



nian H SO = iA SO s^ c ] S - 



Pr,Pt,Vz, £^ VT1 is the Levi-Civita symbol, and Aso = 
6meV [23(1 . The index v = r,t,z labels the spin 
components in the local coordinate system, i.e., S r = 
S x cos cpi + S y sin <pi, S f — S y cos pi — S x sin pi, with pi 
the azimuthal angle of site i (see Fig. [1]) and S x,y,z the 
spin Pauli matrices (with eigenvalues ±1). 

An electric field oriented perpendicular to the tube axis 
affects the electrons in two ways. First, the orbital ener- 
gies are modulated by the electrostatic potential gradi- 
ent. This is described by the on-site energy Hamiltonian 
= eE* Rcos(Lpi)c\^ x Ci^\, where E* is the screened 
electric field, e is the electron charge, and R is the CNT 
radius. This Hamiltonian induces a rearrangement of 
charges on the CNT surface and so, by Coulomb inter- 
action, leads to screening of E. Hence, depends on 
the screened field inside the tube which we find in 
the linear regime to be given by E* = E/j with 7^5, 
in agreement with Refs. 19|-21|. However, the renor- 
malization of the Fermi velocity vp [20[ is found to be 
negligible for the parameters used in this paper. 

Second, an electrostatic potential 0(r) varying on the 
lattice scale induces intra-atomic transitions between or- 
bitals ji and \jl because generally 0(r) \ 7^ 0. Most 
important is the s-p r transition because of two reasons: 
1) It is the only transition directly coupling tt and a 
bands, thus giving rise to a first order effect in the 
s-p r coupling strength. 2) Its strength is determined 

by the unscreened field E and not by E* < E. In- 

(2) 

deed, the induced potential (/>i n d cancels in H K E , i.e. 
(Pr\ </w(0 \s) = 0, as (/>i n d is approximately an even 
function in r about r = R. Based on these arguments, 
we keep only the s-p r transition. The validity of this 
approximation was also verified numerically. The result- 
ing Hamiltonian is H^) = — cE^q cos(^)c| ^Ci s \ + H.c 

where £ = - (Pr\ r \s) = ^ ~ 0.5A with a B the Bohr 
radius and Z ~ 3.2, where we have assumed hydrogenic 
wave functions for the second shell carbon orbitals. 

Effective theory. The microscopic model allows us to 
formulate an effective low-energy theory for the tt band 
near the Dirac points K and K'. As explained, we include 

(2) 

the curvature effects and the s-p r transition H E impor- 
tant for the SOI, and neglect the other inessential interac- 
tions. We have tested this against numerical solutions of 
the full Hamiltonian [Eq. JT}]. We also checked that ad- 
ditional trigonometrically modulated perturbations, such 
as s-pt transitions or sublattice staggered potentials, do 
not change the spectrum qualitatively. 

Hamiltonian ([1]) can be written as H = + H a + 
Hit a, where H n and H a describe the tt and a bands, and 
H ncr the an coupling. For momenta k close to a Dirac 
point 1 1 Hire 1 1 <C 1 1 H n — H a 1 1 . In perturbation theory we 



obtain Hf = H n + H^[H V - H cr ]~ 1 H 7rcr + G((a/R) 2 ). 
Here we keep only terms up to second order in the small 
parameter a/R and the small energies Aso, eE£o which 
must be compared to typical hopping amplitudes ~ eV. 
H ncr ~ + H£q + H E \ where the superscript ttg 

refers to the terms coupling the tt and the a bands. We 
calculate the effective Hamiltonian for the tt band 



rreff rrO , rrcv , rrcv , rrel /o\ 
^tt = H 7r + ^orb + ^SO + **SO> \ 2 ) 

where the last three terms are explicitly listed in Table 
H including numerical values for typical CNTs. Further- 
more, H® = lim/^oo H n = HvpO^qCJi + kra^) is the tt 
band Hamiltonian for flat graphene with periodic bound- 
ary conditions, with r = ±1 labeling the two inequiva- 
lent K and K' points and k the momentum along the 
tube measured from the corresponding Dirac point. For 
semiconducting CNTs, k,Q = (n — rS/3)/R ^ leads 
to a gap 2hvp\kQ\^ where n G Z is the subband index 
and S = (N 1 - AT 2 )mod3 for a (N 1 ,N 2 )-CNT. In the fol- 
lowing, we consider only the lowest subband in metallic 
CNTs defined by k% = Q. 

■^orb d esci> ib es the curvature induced /c-shift of the 
Dirac points [6|, |9|, e.g., K — >• K — Ak cv , with Ak cv = 
(Afc* v , Ak z v ). The shift Afc* v is parallel to the tube and 
can be removed by a gauge transformation shifting the 
origin of k. For non-armchair CNTs, Afc* v ^ and gaps 
are introduced by the curvature H^ h . H^q contains the 
curvature induced SOI @-[||. It contains only S z because 
S r,t depend on cos ip and sin (/?, which average out in the 
if integration. 

On the other hand, H E oc cos<£> which, in 
combination with the SOI terms involving S l = 
S y cos(p — S x sirup, leads to a nonvanishing H^ oc 
S V (J2 Asoe^o / dip cos 2 tp [24]. Since the term propor- 
tional to S r couples only within the a band, it leads to 
negligible higher order corrections. Hence, 



Ht = TeEZS»<T2, 



(3) 



where £ = — £oAso/3Vs P . This is one of our main results. 

Spectrum and helical states. First we focus on arm- 
chair CNTs, assuming that A/c^ v has been gauged away. 
Furthermore, Ak f cv = and f3 = so that the physics 
is completely determined by the interplay of H^ l Q and 
aS z (j\. In Fig. [21 we show the spectrum for a (10,10)- 
CNT in an electric field of lV/nm. For \k\ ^> \a/hvp\, 
H^q aligns the spin in ^/-direction. For the right-moving 
branch (t<t 2 = 1, positive slope) the energy of the S y =t 
state is higher than the energy of the S y =1 state by 
2eE^. For the left-moving branch (r<T2 = —1, negative 
slope) the S y =4- state is higher in energy. Without the 
term aS z <j± the spectrum would be spin-degenerate at 
k = (dashed lines in Fig. [2]). Unlike in usual one- 



dimensional conductors 



these degeneracies can- 



not be lifted by a uniform magnetic field because hy- 
bridization between the crossing bands requires the com- 
bination of spin flip and sublattice hybridization. This 
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TABLE I. The effective Hamiltonian for CNTs. a c± 2.4A is the lattice constant. is the chiral angle (0 
CNTs). <7i,2 are the Pauli matrices in sublattice space. S x,y,z are the spin operators (eigenvalues ±1). (V ss , 
-(6.8, 5.6, 3.0, 5.0, 8.9) eV 1], A SO = 6 meV 23]. The Fermi velocity is v F = VS\Vp P \a/2h ~ 0.95 x 10 6 m/s. 



^ for armchair 



Horb — frVF(Akl v (Ti +tA/Cc V 0"2) * Hvf^^cv = ^Vf 



_ ^pp ( ^pp ^PP ) f ® \ ^ 



^(^pp 



\R 



- cos 3$ 
sin 3$ 



5.4 meV 
R[nm} 2 



- cos 30 
sin 3$ 



Hj& = aS z a! + r/3S z 



V3esAso(Vp P -Vp P j 



lS(V S p) 2 (R/a) 



-0.08 meV -\/3AsoX^cos_3g _ -0.31 me V 

" HVpl + Vg,)(i2/a) ~ i?[nm] 



R[nm] 



cos 30 



flfo = T eEZS y <T2 



S °eo^2xl0- 4 r 



e££ ~ 0.2 meV for E = 1 V/nm 



* j see also Refs. 6-9. 




FIG. 2. (Color online) Energy dispersion e(k) of a (10,10) 
armchair nanotube (implying R — 0.67nm). The solid (red) 
lines show the analytical results [see Eq. J4} and Tab. U and 
the dots show numerical results obtained from H in Eq. JT]). 
The axial /c-shift Ak z v has been removed in both, the numer- 
ical and the analytical spectrum. The arrows correspond to 
the S y projections and are reversed for E — » —E. The field 
strength is E — 1^ so that the splitting 2eE£ > c± 0.4meV 
and the gap 2\a\ ~ 0.24meV. The dashed lines indicate the 
spectrum for the case a = with the spin-degeneracies at 
k = 0. The dashed gray (light blue) lines indicate chemical 
potentials at which only helical modes exist [see Eq. J5J)]. 

is, however, caused precisely by aS z ai, which is gener- 
ated by virtual transitions to the a band that result in 
the simultaneous spin and sublattice hybridization. As 
a result, a gap of size 2\a\ is opened at each degeneracy 
point. The resulting spectrum, shown in Fig. [2j at the 
K point has four branches, the subbands |ra, fc), given by 



e(k) = ±eE£ ± V 'a 2 + (hv F kf 



(4) 



An equivalent spectrum exists at K'. The spin orienta- 
tions on the branches for \k\ ^> \a/hvF\ are identical at 
both Dirac points (arrows in Fig. [2]). For general fc, the 
S y expectation value in state \m,k) is given by 



(m, k\ S y |m, k) = ±k/ v / (a/hv F ) 2 + k 2 , 



(5) 



where for eE£ > the + corresponds to subbands m = 
1,4 in Fig. [2j and the — to m = 2,3 (and vice versa 
for eE£ < 0). Note also that the expectation values of 



S x and S z in all states give zero, so that only (S y ) ^ 
0. In this sense, the states are always perfectly spin- 
polarized, even though the measured total spin is not 
unity. The bands crossing the chemical potentials /ii 
and fi2 indicated in Fig. [2] have (S y ) ~ ±0.95. We also 
note that electron-electron interactions generally lead to 
an enhancement of the gap 2\a\ [141 ]. 

Fig. [2] shows the analytical spectrum Eq. (|4j) for an 
armchair CNT in comparison with a numerical diagonal- 
ization of Hamiltonian (pQ). The qualitative features of 
the spectrum are well preserved by the effective theory. 

If, in an armchair CNT, the chemical potential is tuned 
to iii or fi2 (see Fig. [2]), the remaining conducting modes 
are helical, i.e., the direction of motion is coupled to the 
spin direction. In the present case, the spin points along 
E x v, where v = ±^pz for right and left movers, re- 
spectively. In particular, this implies that E — >> — E 
also reverses the helicity, thus inverting the spin filter- 
ing. We note that the helical modes are stable against 
small deviations from the (TV, 7V)-CNT (armchair) case 
with chiral angle 9 = ^. The additional terms j3S z 
and hvp/Sk^ai, which appear for 9 ^ ^, partially align 
the spin in z-direction and open gaps at the zero-energy 
crossing points. We find that for metallic chiral CNTs, 
e.g. with (N + 3, N) and N ~ 10 - 20, that are close to 
the armchair limit, good spin polarization ((S y ) ~ 90% 
and (S z ) < 20%) can still be obtained (see also Fig. [3]). 

Valley suppression. In chiral (N + 3/,iV)-CNTs, with 
I = 1,2, it is possible to mostly restore the armchair 
spectrum and spinor properties for one Dirac point by 
the further application of a magnetic field B z along the 
tube. As mentioned above, 9 < ~ results in cos(3#) ^ 0, 
thus leading to two additional terms in the Hamiltonian: 
a transverse /c-shift Hvp^k^ai and an effective Zeeman 
field rf3S z . These terms have opposite signs at differ- 
ent Dirac points. The field B z leads to terms of the 
same form, yet with equal signs at both Dirac points, 
so that the chirality-induced cos(3#)-terms can be can- 
celed at one of the Dirac points, whereas they are dou- 
bled at the other. Indeed, the orbital effect of B z adds 
Akg = 7rB z R/&o to Afc* v , with $0 the magnetic flux 
quantum. The Zeeman effect of B z adds hbB z S z to 
r/3S z . Due to the different radius (R) dependencies of 
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FIG. 3. (Color online) Chiral (17,14)-CNT with electric field 
E = lV/nm (eE£ ~ 0.2meV). (a) The spectrum and (b) the 
spin expectation values at the K/K' point for the |2, k) sub- 
band and magnetic field B z = 0. For B z = 0.8 IT, the bands 
(c) at K' (dashed lines) are gapped, while the spectrum at K 
(solid lines) has the same form as in the armchair case [Eq. 
dU)]. The size of the gap is 2\a\ = 0.16meV. The spin expec- 
tation values (d) at K for |2, k) follow closely the armchair 
case [see Eq. ©]. 



the Zeeman and orbital terms, R and B z can be cho- 
sen such that both cos(3#) terms in Table U cancel at 
one Dirac point, provided that Aso > [25] (see Fig. 
[3]). However, since R cannot be chosen continuously, the 
cancellation is perfect only for one of the two terms. The 
Zeeman term can be removed at K with B z = —/jlb/ 'ft, 
but small gaps will remain at energy s = 0. On the other 
hand, if the e = gaps are to be closed, B z must be 
tuned such that Ak f B + Afc* v = 0. The small residual 
Zeeman term /3* = f3 + [i B B z (|/3*| <C \ ft\) then leads to 



e = ±V(/?*) 2 + (e£0 2 ± Vi^Fk) 2 



a 



(6) 

±/3*/eE£ 



and to a small spin-polarization (S z ) 
((S x ) = in all cases). 

An illustrative example is the (17,14)-CNT with B z c± 
0.81T, for which the orbital and Zeeman cancellations 
work particularly well. At K the spectrum and the spinor 
properties (m, k \ S XjVjZ \m, k) of an armchair CNT are re- 
stored, while at K' the curvature-induced gap hvpAk 1 ^ 
is amplified by a factor of 2 (see Fig. [3]). This ampli- 
fication is sufficient to remove all states of K' from the 
relevant energy range so that only K contributes a single 
pair of helical modes at the chemical potential 

We note that as an immediate consequence of the SOI 
induced gaps the conductance of the CNT is reduced by 
a factor of two, and by an additional factor of two if 
the valley degeneracy is lifted. As mentioned, helical 
modes can be used as spin filters, Cooper pair splitters, 
and allow for Major ana fermions at the CNT edges if the 
latter is brought in contact with a superconductor. These 



properties, together with the all-electric control, make 
CNTs attractive candidates for spintronic and quantum 
computing applications. 
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